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Abstract 

We investigate the presence of topological twistons in crystalline polyethylene. We describe 
crystalline polyethylene with a model that couples the torsional and longitudinal degrees of free- 
dom of the polymeric chain by means of a system of two real scalar fields. This model supports 
topological twistons, which are described by exact and stable topological solutions that appear 
when the interaction between torsional and longitudinal fields is polynomial, containing up to the 
sixth power in the fields. We calculate the energy of the topological twiston, and the result is in 
very good agreement with the value obtained via molecular simulation. 
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The existence of twistons in crystalline polyethylene (PE) was postulated two decades ago, 
and refers to a twist of 180*^ that extends smoothly over several CH2 groups in crystalline PE, in 
the plane orthogonal to the chain direction, with the corresponding CH2 unit length contraction 
along the polymeric chain. These twiston configurations appear in crystalline PE as a result of 
its large torsional flexibility, and may contribute to elucidate some of its properties, in particular 
the dielectric a relaxation 0, |^, |^, |, |^, 



There are some interesting models of twistons in crystalline PE p, ^, |T^, [11], [T2|. The works 
1^ are almost simultaneous to the work |jl3|, which introduces solitons to describe conductivity 
in polyacetylene via distortions of the single-double bond alternations. In the PE chain, however, 
the bounds are always single bounds, which require at least one bosonic degrees of freedom 
to describe the torsional flexibility of this unsaturated polymer. Despite these two decades of 
investigations, we believe that the issue of topological twistons playing some role in explaining 
properties of the crystalline PE chain is still incomplete, requiring further investigations both 
in the theoretical and experimental grounds. This is the main motivation of the present work, 
where we follow an alternate route to topological twistons in PE to bring new facts to the former 
theoretical investigations. This new route was introduced in Ref. [l^, and here we complete the 
investigation, including the calculation of the energies of the exact topological twistons. 

We start our investigations by first reviewing the basic features of the several distinct mechan- 
ical models used to describe twistons in crystalline PE. The most important ones are described 



in Refs. P, |T0|, |TT|, . In the pioneer work p[ the author considers a system which couples the 
torsional and translational degrees of freedom. In Ref. P] the authors consider a simpler system, 
describing only the torsional motion along the crystalline chain, and this is also considered in 



the subsequent work [|10l . It is only more recently ITT], that one includes interactions between 
radial, torsional and longitudinal degrees of freedom. In this case one uses cilindrical coordinates 
to describe a generic CH2 unit via (r„, z„), which correspond to the three degrees of freedom 
of the rigid molecular group. A simplification can be introduced, and concerns freezing the r„'s, 
so that the radial motion is neglected. In |ll[] one further ignores the translational degree of free- 
dom, the Zn coordinates, to get to a simple model described via the torsional variable that in the 
continuum limit can be taken as 9n{t) — > 6{z,t). The model reproduces the double sine-Gordon 
model, according to the assumptions there considered to describe the intermolecular interaction. 
The other more recent work |]T^ on twiston in crystalline PE gives another step toward a more 



realistic model. This is the first time the radial, torsional and longitudinal degrees of freedom 
are simultaneously considered to model twiston in crystalline PE. The model is very interesting, 
although it is hard to find exact solutions and investigate the corresponding issues of stability. 
The problem engenders several intrinsic difficulties, which have inspired us to search for an alter- 
nate model, in the form of two coupled fields belonging to the class of systems investigated in the 
recent works 0, [T^, [!§ . 

The basic assumptions introduced in the former models for twistons in crystalline PE may be 
described considering cilindrical coordinates. The Lagrangian presents the usual form L = T — U, 
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where 



n 



U — t/intra + fAnter (2) 

Here m is the mass associated to the molecular group CH2, and ?7intra and U\^ter are potentials 
used to model the intramolecular and intermolecular interactions in the crystalline environment, 
respectively. The intramolecular potential can be considered as 

f/intra = ^ E ^1 (^"+1 " <^nf + j E ^2 " Z^f + • • • (3) 

^ n ^ n 

where Ki and K2 are coefficients related to the harmonic approximation for torsional and longitu- 
dinal motions, respectively. The intramolecular potential may contain derivative coupling between 



the torsional and longitudinal motions. In this case we should add to Uintra the contribution []TT| 



i E ^3 {dn+l - 0nf " ^n) (4) 

^ n 

However, instead of the above coupling we can consider derivative coupling between the radial 



and longitudinal motions. In this other situation we add to Uintra the contribution |12 



1 1 

9 E ^4 (r„ - ro)^ + - E ^5 {r-n+l - r-n) {Zn+1 " ^n) (5) 

The above terms (§) and (|]) are two among several other possibilities of introducing derivative 
coupling between the torsional, longitudinal and radial degrees of freedom. We shall not consider 
such possibilities in the present work, although in one shows a route for taking derivative 



coupling into account. As we are going to show below, we consider the standard harmonic 
approximation for Uintra in order to follow the basic steps of the ffist works |10| on twistons 
in crystalline PE. 

The second potential in (0) is U inter- It responds for the intermolecular interactions and is 
usually given in the form f/inter = Hn \Uo{On) + f^i(^n) Ui{zn)\- Here Uo{9n) and f/i(6'„) are used 
to model torsional mobility and Ui{zn) to describe the longitudinal motion along the chain. In 
the works |^, |TT|, after freezing radial and translational motion, the above intermolecular 
potential is described by the Uo{9n) contributions. We can get to models for the torsional motion 
alone, and in the continuum limit they may be described via the sine-Gordon potential 10 



Ai [1 — cos(2 6')], or the polynomial potential 6*^ + Bi9^, or yet the double sine-Gordon 
potential |Tl[| A3 [1 — cos(2 6')] + i?2 [1 — cos(4 6')]. Here and Bi are real constants, used to 
parametrize the corresponding interactions. Evidently, the above potentials lead to different 
models for the torsional field, and are introduced to account for the specific motivations presented 
in the respective works |^, |TD|, |TT]. In the more recent work |T^, one considers coupling 
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between the radial, torsional and longitudinal degrees of freedom, but the analytical solutions 
there obtained are found under assumptions that ultimately decouple the system. 

The above models show that the basic idea introduced in Ref. [Q has survived along the 
years, although there have been interesting quantitative contributions to investigate the presence 
of twistons in crystalline PE. In particular, in Ref. [|l^ one includes the most relevant degrees 



of freedom when one considers the CH2 group in the form of rigid molecular group along the 



crystalline chain in crystalline PE. However, in the model considered in [|I^ we could not fully 
understand the reason for not considering harmonic interactions between neighbor radial coordi- 
nates, while taking into account interactions between radial and longitudinal degrees of freedom 
in the intramolecular potential. For this reason, we think we can introduce another mechanical 
model for the polymeric chain, where we modify some assumptions presented in the former works 



11] , |T2| . The difficulties inherent to the problem of describing topological twistons in crystalline 
PE bring motivations for simplifying former assumptions, with the aim of offering an alternate 
model that presents exact solutions for twistons in crystalline PE. Toward this goal, let us use 
cilindrical coordinates to describe the molecular groups under the assumption of rigidity. We 
start with the kinetic energy ([^), rewriting it in the form 

T=lmrlY.[^l+{0 Xl + P'r}j (6) 

Here we have set 0„ = 6'„ - [1 - (-l)"](7r/2), Xn = {zn - nc)/c and p„ = (r„ - ro)/ro, where 
ro is the equilibrium position of the radial coordinate and c is the longitudinal distance between 
consecutive molecular groups. Now 0„, Xn and p„ are all dimensionless variables, and in the 
continuum limit they can be seen as real fields (j){z,t), xi^^t) and p{z,t). Before going to the 
continuum version of the PE chain, however, let us reconsider the intramolecular potential given 
by Eq. (^. We use the harmonic approximation to write 

111 

f^intra = 7: ^ h (0„+i - + 77^2^1 (^ri+l " Xnf + 77 J2 {Pn+1 " Pnf (7) 

^ n ^ n ^ n 

where kt, ki and k^. are spring-like constants, related to the torsional, longitudinal and radial 
degrees of freedom, respectively. 

The harmonic interactions present in the intramolecular term (|^) makes the dynamics to 
appear as the dynamics of relativistic fields, in the same way it happens with the standard 
harmonic chain. We use © and (0) to write the following Lagrangian density for the continuum 
version of the mechanical model for crystalline PE 

Im 2 ( dcpy 1 ( d(i)V Im 2 / dxV 

^'^(S) -^--('^'^'Z^) 

The quantity m/c identifies the mass density along the chain, and ktc = Kt, hic = 

are Young parameters related to the torsional, longitudinal and radial motion, respectively. 
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The above mechanical model is still incomplete, but it contains the basic assumption that we 
are dealing with an harmonic chain, and deviation from the harmonic behavior is to be included 
in Vinter- Although in this case we can not introduce any other derivative coupling, we still have 
the freedom to specify Vmter and so introduce nonlinearity via the presence of the surrounding 
environment in the crystalline material. This is the model we keep in mind to introduce the 
following field theoretic considerations. 

We follow the lines of the former mechanical model, which lead us to introduce a field theoretic 
model that contains three real scalar fields. The Lagrangian density describing the fields (p = 
(f){x,t), X = xi^yt) and p = p{x,t) was introduced in Ref. Jl^. The model is defined by the 
potential, which is supposed to have the form 

v{<P,x,p) = IhI + ^h'^ + ^hI (9) 

Here Hip = dH/d(j) and so forth. H = H{(j), x, p) is a smooth but otherwise arbitrary function of 



the fields. This restriction is introduced along the lines of the former investigations |15, |T^, 17], 
and leads to interesting properties, such as the ones explored below. 

We focus attention on the crystalline PE chain. In this case it is a good approximation 
P, |10|, [ll], |12| to descard radial motion in the PE chain. This simplification leads to a system 



of two fields, describing torsional and longitudinal motions simultaneously. However, we first 
consider the simpler system, described by the torsional field alone. In this case, in accordance 
with the Refs. |Tl| , |T2[, investigations on molecular simulation allows introducing the following 
torsional potential 

V,{<P) = \\'<P\<P'-'n^f (10) 

Fortunately, this potential is generated by the function Hi{(j)) = (1/2) (0^/2 — tt^). Also, it 
has three degenerate minima, one at = and the other two at 0^ = vr^. 

We use the potential Vi (0) to get the masses of the elementary excitations around = and 
= ±7r in the form m0(O) = |A| tt^ and m^{±Tc) = 2 |A| tt^. These results identify an asymmetry 
in the spectra of excitations of the torsional motion around the minima = and 02 = vr^. This 
asymmetry appears in consequence of the polynomial potential (p!0|), and is small for small A. It 
is related to the asymmetry between the well at = 0, and the well at 0^ = tt^. Since the mass 
of the field corresponds to the minimum energy necessary to excitate elementary mesons into the 
system, we realize that the value |A| vr^, the difference 2 |A| vr^ — |A| tt^ may be seen as the energy 
for the field go from = to = ivr, that is the energy to overcome the torsional barrier in 
this simplified model. 

To get to a more realistic model we couple the torsional field to the longitudinal motion along 
the chain. We model the presence of interactions by extending the former function i/i(0) to 
^^2(0,x) given by 

i/2(0, X) = ^ A 02 (^02 - n') + ^p4>'x' (11) 
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This gives the system 



where 



V.{^^ X) = - A' + A/.</>^(</>^ - ^')x' + + ^/i W (13) 



We are using natural units, as in Ref. |14|. The above potential presents interesting features. 
For instance, V2(0, 0) = Vi((/)), which reproduces the torsional model ^^(0) when one freezes the 
longitudinal motion. Also, V2(0,x) = and 

^2(±vr,x) = ^/iW + ^/iW (14) 

We can evaluate the quantity S^V j d(\)dx to see that it contributes with vanishing values at the 
minima (0,0) and (ivr, 0). This shows that the spectra of excitations of the torsional motion 
around the ground states are unaffected by the presence of the longitudinal motion. Thus, we 
can use Vi(0) to investigate the behavior of the torsional motion around the equilibrium config- 
urations. 

The masses of the field are now given by m0(O, 0) = |A|7r^ and m(^(±7r, 0) = 2|A|7r^, around 
the minima (0, 0) and (±7r, 0), respectively. Accordingly, for the x fi^ld we see that it is massless 
at (0 = 0, X = 0), and at (±7r, 0) the mass is m^(±7r, 0) = |/i|7r^. These results identify an 
asymmetry in the spectra of excitations of both the torsional and longitudinal motion around the 
minima (0,0) and (7r^,0). This asymmetry appears in consequence of the polynomial potential 
(p!3D, and is small for small parameters A and [i. These results allow introducing the ratio m^jm^ 
between the masses of the and x fields - see Ref. fl^ . 



The topological solutions connect distinct, adjacent minima of the potential. The energy 
corresponding to the classical configurations can be written in the general form JIT 



Eij = \H{^i,Xi,pi) - H{^j,Xj,Pj)\ (15) 

where {4>i,XiyPi) and {jPj^XjiPj) stand for two vacuum states, that is, two adjacent points i and 
j in the field space (0, x, p) that minimize the potential. 

Let us first consider the case of a single field, the field that describes torsional motion along 
the polymeric chain. We use former results to write the equation of motion for static configuration 
in the form 

^ = AW-^^)(30^-vr^) (16) 
This equation is solved by solutions of the first-order equation 

f = A0(0^ - TT^) (17) 

az 
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There are topological twistons, given by [|T5] 



cf>\^\z) = ± TT Y/(l/2)[l-tanh(A7r2^)] (1^ 



Here we are taking 2; = as the center of the soliton, but this is unimportant because the 
continuum, infinity chain presents translational invariance. The sign of A identifies kink and 
antikink solutions, connecting the minima and vr or and —it. These solutions are stable and 
can be boosted to their time- dependent form by just changing z to = (z — vt) / (1 — v'^Y^'^ . This 



model can be seem as an alternate model to the ones introduced in the former works [§, ^, [ly, |Il 

The amplitude of the torsional field is tt, which is the angle the chain rotates to form the 
twiston. The width of the twiston, L(i), which is the length along the chain where the angular 
position of CH2 groups appreciately deviates from the crystalographic positions, is inversely 
proportional to the quantity jAlvr^. We can also get the energy corresponding to the static 
twiston. We use Eq. (|15D to get the value 

E^t) = \\\\7i' (19) 

We now consider the model that describes interactions between the torsional and longitudinal 
fields. The equations of motion for static fields (p = (f){z) and x = x(^) ^I's given by 



fx 

dz^ 



AV(0' - vr^)(302 - vr^) + 2A^0(202 _ vr^)^^ + + l)x' (20) 

2A/i02(02 - vr2)x + 2/i202^3 ^ ^2^2^ (21) 



Although there is no general way of solving these equations, we recognize that they follow from 
the potential in Eq. (0), defined via the function introduced in Eq. ([TI|), and so they are solved 
by 

^ = A0(02 - vr^) + (22) 
az 

^ = (23) 
dz 

which are first-order differential equations, easier to investigate. 

To find explicit solutions we use the trial orbit method introduced in Ref. IT^]. We consider 
the orbit 

A(02-7r2)+/ix' = /i(0'-vr2) (24) 



We note that this orbit is compatible with the first-order Eqs. ( P2|) and (^). Also, from Eq. (p^) 
we get 

(t^fl (z) = ± TT V'(l/2)[l-tanh(^7r2^)] (25) 
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This result and the orbit (plf) are now used to obtain 



X[g(^) = ±vry^-1 V(l/2)[l + tanh(/ivr2^)] (26) 
These solutions are valid for A//i > 1 and are similar to the solutions found in Ref. [|12| to 



describe the torsional and longitudinal degrees of freedom that describe topological twistons in 
the crystalline PE chain. 

The amplitude of the twiston is still vr, while the amplitude of the longitudinal motion is given 
by 7r[(A//i)-l]^/2_ This result requires that A/yU > 1, which is compatible with the investigation of 
Ref. [|T^. In this more sofisticated model the width L(^t,i) of the topological twiston is proportional 
to l/(|yu|7r^). It depends inversely on /x. We compare L(t) and L(^t,i) to see that L(^t,i) > -^(t) since 
A//i > 1 for the topological twiston of the model of two coupled fields. This result is new, 
and shows that the presence of the longitudinal motions contributes to enlarge the width of the 
topological twiston. 

Another result follows after calculating the energy of these solutions. We use Eq. (^) to get 
E(^t,i) = -£'(4) = (l/4)|A|7r'*, which equals the value of the energy of the simpler model, where one 
discards the motion of the longitudinal field. This result shows that although the more general 
model changes some of the features of the simpler model, which describes only the twiston field, 
it does not change the energy of the twiston. We understand this result as follows: the first-order 
equations (|^) and (|23D also present the pair of solutions 



^fj^ {z) = ±n ^(l/2)[l-tanh(A7r2^)] X(t,0 (^) = (27) 



This pair of solutions and the former one, given by Eqs. (^) and (pq), are at the same topological 
sector and present the very same energy, given in Eq. (??). However, when one sets x — >■ in 
the coupled model, the system changes to the simpler model, and so the energy of the pair (P^) 
is necessarily equal to the energy of the twiston in the single field system. This fact explains our 
results, and shows that the torsional energy is the main quantity to calculate the energy of the 
topological twiston. We use this point of view to rewrite the energy as = (1/4) (|A|7r^) vr^. 
We have already identified |A|7r^ and 2|A|7r^ as the masses of the twiston field, which show that 
when the field varies from to ivr, that is when a twiston is formed, one changes from the 
energy |A|7r^ to the energy 2|A|7r^, and this requires the value |A|7r^. We then identify this value 
with the energy for twiston formation along the polymeric chain. According to Ref. 0] , the energy 
contribution of the localized twisted region to the criation of the twist defect is 7.3 Kcal/mol. In 
fact, in Fig. [6] and Table II of Ref. we see that Uq = 9.8 — 2.5 = 7.3 Kcal/mol, which is to 
be regarded as the contribution of the localized twisted region to the criation of the twist defect 



20|. We then change |A|7r^ —>■ 7.3 in the energy to obtain 

= 17.99 Kcal/mol (28) 

This is the energy of the topological twiston, and is in good agreement with the energy values 
of 18.01 Kcal/mol 0, 18-19 Kcal/mol |Tg, and 17.2 Kcal/mol fT^, obtained using different 
numerical simulations and models. 
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We conclude this letter recalling that we have investigated a system of two coupled real scalar 
fields to model topological twistons in crystalline PE. This model describes no radial motion, but 
it couples the torsional and longitudinal degrees of freedom in a very interesting way. We have 
found exact solutions, which engender several features, and here we offer the following remarks. 
The limit — > A transforms the solutions (^51) and (^) into the solutions (0) of the former case, 
that describes the torsional motion alone. The solutions of the model of a single field present 
width proportional to 1/|A|, and for the two field model it is proportional to l/|yu|. The width of 
the solutions of the two fields is exactly the same, in agreement with the topological features of 
the solutions, and with the orbit (p^ , used to solve the coupled equations (p2D and (|^). This 



result is intuitive, since one expects that when the torsional motion completes the 180° rotation 
and returns to its crystalographic position the longitudinal motion should simultaneously return 
to its crystal register. The amplitude of the torsional field is given by the solution (|25|) and 
is TT, in agreement with the model we use for the twiston configuration. The amplitude of the 
longitudinal motion is given by the solution (^) and is vr [(A//i) — 1]^/^. In the PE chain we have 
to set this to unit, to make it compatible with c, the full longitudinal motion. This picture follows 
in accordance with the fact that crystalline PE presents degenerate ground states, obtained from 
each other by a rotation of 180° or by a translation of c along the polymer chain. 

We have also obtained the energy of the twistons. It is = (1/4) |A| vr^. We have used |A| vr^ 
to identify the mass difference for the torsional field in the minima = vr and = 0. This and 
results of Ref. [|l| allow getting |A| tt^ = 7.3 Kcal/mol, which gives the energy of the topological 
twiston as 17.99 Kcal/mol, in good agreement with values known in the literature. 



The results presented in this work completes the former investigation [Q. They show that 
the approach of using systems of coupled fields and the corresponding field theoretic analysis to 
describe topologically non-trivial excitations in continuum versions of polymeric chains seems to 
work correctly. The procedure describes interesting aspects of the problem, and allows obtaining 
the energy of the topological excitation in a direct way. We believe that similar polymeric chains 
can also be investigated by similar systems, and this makes us to think on modelling topological 
twistons for instance in the family of systems where one changes some CH2 groups by oxygens 
periodically, to make chains with the basic units CH2 — O, CH2 — CH2 — O, CH2 — CH2 — CH2 — O, 
etc. Despite the presence of oxygen the bounds are still sigma bounds, and the torsional motion 
seems to be similar to the PE chain. Thus, we may use a twiston model to explore properties 
of the family {CH2)n — O, in particular in the case of CH2 — CH2 — O, the Poly(oxyethylene), 
POE. This and other related investigations are presently under consideration. 
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